Journal of Theoretical Biology

This paper examines the impact of Cancer Care Ontario's regionalization of lung cancer surgery on blood transfusion use for thoracic surgery at The Ottawa Hospital. Using aggregate data scraped from the paper ''Changing trends in blood transfusion: an analysis of 244,013 hospitalizations'' by Shehata et al. (2014) we estimated the effect regionalization on the log mean number of blood transfusions among thoracic surgeons (who treat lung cancer), utilizing other surgical specialties as controls. We found that the average number of blood transfusions increased significantly among thoracic surgeons following the opening of a new cancer care centre at the hospital in December 2009, while control specialties showed no parallel change. Contextualizing these results in a modified ''population level'' treatment threshold model, it is argued that the main driver of this increase in blood transfusions among thoracic surgeons was due to the influx of less skilled thoracic surgeons to The Ottawa Hospital, and the closing of thoracic surgery units at smaller hospitals which was part of the broader regionalization process.

The -Bonacci weighted sequence space {ell}(,) is introduced for [&ge;] 2, 1 [&le;] < {infty}, and > 0). It is proved to be a Banach space which is isometrically isomorphic to {ell}, with an unconditional Schauder basis. The space is also reflexive and uniformly convex. Its strict inclusion chain, continuous dual, and asymptotic truncation estimate of order - / are also established. Based on this framework, Scaled &#119898-BonacciNet is proposed as a lightweight CNN that rescales Fibonacci-Net channels using the dominant &#119898-Bonacci root and employs depth wise separable convolutions. Tribonacci-Net achieved 96.09% accuracy and a 0.9581 macro F1-score on the 15-class Brain MRI dataset using 0.85 million parameters and 0.366 GFLOPs. Pentanacci-Net attained 95.87% accuracy on the 44-class dataset while reducing parameters and FLOPs by 47.10% and 44.48%, respectively. On the CKD ultrasound dataset, Hexanacci-Net achieved 83.75% accuracy and a 0.7941 macro F1-score, with corresponding reductions of 50.37% and 45.10%. These findings support dominant-root scaling as a structured approach to lightweight CNN design.

We develop an age- and immunity-structured COVID-19 model with setting-specific contacts, calibrated to Ontario incidence and hospital admissions before Omicron. We simulate counterfactual vaccination strategies and contact restrictions in schools, workplaces, and other community settings, evaluating projected severe-case demand against an illustrative healthcare-capacity benchmark. Vaccination combined with workplace and community contact restrictions keeps simulated demand within capacity. In contrast, relaxing restrictions in either setting produces large surges and repeated capacity threats, even when school restrictions remain. School-only changes produce smaller effects on severe disease. These findings highlight the value of vaccination and setting-specific contact restrictions that prioritize workplace and community settings for pathogens with adult-concentrated severe outcomes.

The aim of this study was to develop a framework for assessing the incremental value of decisions on drug amount, timing, input structure, cessation and restart. The LOGOS (Longitudinal Optimization of Global Outcomes and Strategies) approach combines matched nested regimen classes, mechanism-specific benchmarks and evidence requirements in a single longitudinal comparison. Thus, the contribution of additional regimen freedom is assessed separately from improvements due to different resources, constraints or evaluation horizons. Nested feasible classes have nondecreasing objective suprema, whereas a strict gain requires a feasible witness above a valid bound for the smaller class. The approach is illustrated through analytical examples of periodic smoothing, matched desensitization, relaxed exposure allocation, finite-course loading and safety after cessation. In a normalized finite Emax model, integrated benefit increased from 0.629606 for three equal boluses to 0.672945 after optimization of amounts, 0.677043 for an amount-and-time candidate and 0.679553 for a globally certified loading-maintenance-coast solution. The adjacent observed gains were 6.88%, 0.61% and 0.37%, respectively, and all schedules satisfied the full-tail injury bound. In addition, a periodic factorial example separated response curvature from a joint recovery and repair process, including a null case and results dependent on feasibility constraints. The LOGOS framework allows assessment of whether additional regimen flexibility improves benefit, provides no added value or restores feasibility. These dimensionless examples support conditional mathematical and model-based conclusions; clinical superiority and safe restart in patients have not been established.

Brazil experienced its worst dengue epidemic in 2024, with an increase in severe outcomes among elderly individuals, highlighting the need for vaccination strategies informed by population immunity. We reconstructed state- and age-specific immune profiles by fitting an age-structured catalytic model to dengue surveillance data from 2000-2025. This inferred immune profile was incorporated into an SIR-type transmission model to project disease burden and to evaluate the public health impact of vaccine allocation strategies. We observed substantial spatial heterogeneity in accumulated immunity, which was closely associated with projected burden and optimal allocation. Strategies prioritizing states with lower historical incidence and larger susceptible populations yielded greater benefits. Under limited vaccine supply, targeting children resulted in the greatest reductions in DALYs, whereas targeting elderly populations achieved smaller DALY reductions but larger reductions in mortality. We further showed that accelerated vaccine supply markedly increased impact, while allocation choices remained critical when availability was constrained.

Estimating the correlation between subject-specific rates of change is a common objective in longitudinal studies. A natural approach is to fit ordinary least squares regressions per subject and correlate the resulting slope estimates directly. Under complete data, this slopes-first estimator is algebraically equivalent to a contrast-based plug-in estimator derived from the joint variance-covariance matrix of the two signals across time. Under independent dropout, however, the slopes-first approach implicitly applies heterogeneous weights across missing-data patterns, introducing bias whenever the cross-covariance between the two signals is not constant over time. We formalise this result and propose an unbiased alternative that estimates the joint covariance matrix before applying the contrast. We further extend our results to include covariates as well. Simulations confirm that the slopes-first estimator is biased under non-constant cross-covariance while the proposed alternative remains unbiased. Further, we show in a real-world example that the two estimators can lead to different scientific conclusions and recommend that researchers working with longitudinal data inspect the missing data pattern to make an informed decision about the estimator of choice.

Background: The C-reactive protein-albumin-lymphocyte (CALLY) index integrates systemic inflammation, nutritional status, and immune function. We evaluated its association with survival in primary lung malignancies. Methods: Embase, PubMed, Scopus, and Web of Science were searched from inception to 18 September 2026. Risk of bias was assessed using the Quality In Prognosis Studies tool. Hazard ratios (HRs) were synthesized using random-effects models with restricted maximum likelihood estimation and Hartung-Knapp adjustment. Sensitivity analyses, subgroup analyses, and meta-regression examined robustness and heterogeneity. Results: Eleven studies contributed 13 cohorts to the overall survival (OS) analysis and 8 cohorts from seven studies to a combined recurrence/progression-free survival (RFS/PFS) analysis. Lower CALLY was associated with poorer OS (HR 1.78, 95% CI 1.49-2.12; I2=71.6%) and poorer RFS/PFS (HR 1.60, 95% CI 1.40-1.84; I2=0%). Associations persisted in leave-one-out analyses. No statistically significant subgroup differences were detected. Funnel-plot asymmetry suggested small-study effects for OS and RFS/PFS. Seven studies had moderate and four had high overall risk of bias. The certainty of evidence was rated as low for both OS and RFS/PFS, primarily due to risk of bias, heterogeneity, and endpoint indirectness. Conclusion: Lower CALLY was associated with adverse survival outcomes in the current synthesis. Methodological limitations and variation in biomarker measurement and endpoint definitions restrict interpretation. Prospective validation using standardized measurements and clinically coherent outcomes is needed to establish incremental prognostic value and clinical utility. Keywords: CALLY index; lung cancer; prognosis; overall survival; systematic review; meta-analysis

Vegetation plays a critical role in shaping fine-scale tick density by structuring local microclimates, habitat resources, and host interactions that largely determine tick survival and activity, with implications for tick-borne disease transmission. Tick burden in the United States continues to increase, with annual reported tick-borne disease cases more than doubling from 2004 through 2019. These trends call for a broader scope of tick-mitigation strategies that go beyond traditional prevention measures that focus on personal protection; and there is an opportunity to inform habitat-level interventions. Here, we quantify fine-scale vegetation structure across 48 subplots at tick surveillance sites in Western Massachusetts to assess local habitat features associated with Ixodes scapularis density and potential for human-tick contact risk. In assessing understory and overstory vegetation attributes, this study identified life stage-specific differences in predicted tick density, with higher adult densities in areas with deeper leaf litter and presence of invasive plant species and nymphs associated with high presence of ferns. Both adult and nymphal Ix. scapularis were predicted at higher densities in relatively low basal area forests (18.9 to 26.9 m2/hectare). This study design provides a simple, low-cost framework for organizations looking to incorporate tick-informed forest measurement ecology into their public health tick surveillance protocols with minimal training.

Background: Conventional performance metrics do not indicate whether artificial intelligence (AI)-platforms and radiologists are right or wrong on the same examinations, nor whether their non-agreement is more often beneficial or harmful. Objective: To develop, illustrate, and apply directional metrics quantifying non-agreement between radiologists and AI-platforms against a common reference standard, and to adapt a fragility-based framework for quantifying the statistical robustness of these metrics. Methods: We defined Positive Non-Agreement (PNA, AI-platform correct, radiologist incorrect), Negative Non-Agreement (NNA, radiologist correct, AI-platform incorrect), and Net Non-Agreement (NENA, their difference). Metrics were illustrated in three hypothetical scenarios of 1,000 examinations each and applied to three previously published AI-radiologist studies on pulmonary embolism and mammography. Bootstrap confidence intervals were calculated, and a reclassification fragility index (FI) quantified the minimum case reclassification needed to remove statistical support for NENA's direction. Results: In illustration, NENA was positive in two scenarios and negative in one. In the three published studies, NENA ranged from -4.18% to +10.74% and decreased consistently with increasing radiologist experience (junior residents, +6.27%; consultants, +1.45%), a gradient driven primarily by declining PNA with modest increase in NNA. The reclassification FI was defined for all three illustrative scenarios (range, 0.10%-1.70%) and for six of the seven comparisons from these studies (range, 0.34%-4.31%), but not for the weakest comparison, whose confidence interval for NENA included zero. Conclusions: PNA, NNA, and NENA provide a directional, case-level complement to conventional performance metrics for AI-radiologist non-agreement. The accompanying reclassification FI offers a statistically grounded, generalizable method for quantifying their robustness.

Financial well-being (FWB) encompasses objective and subjective financial status and may influence health. We developed and validated a parsimonious FWB index using data from 23,006 Health and Retirement Study (HRS) participants aged 50 years and older from 2006 to 2022. The FWB-5 Index was derived using correlated two-factor confirmatory factor analysis accounting for household clustering. Indicators included household wealth, household income-to-poverty ratio, monthly financial strain, satisfaction with one's financial situation, and perceived financial mastery. We assessed convergent and discriminant validity using the Consumer Financial Protection Bureau Financial Well-Being (CFPB-FWB) Scale and Satisfaction with Life Scale, respectively. Associations with pre-baseline socioeconomic factors and subsequent depressive symptoms and health conditions were examined to assess construct and predictive validity. Model fit indices supported the two-factor structure (robust comparative fit index = 0.99; robust root mean square error of approximation = 0.05). Correlations were stronger with the CFPB-FWB Scale (Spearman P = 0.80) than with life satisfaction (Spearman P = 0.43). Associations with socioeconomic factors and subsequent health outcomes were in the hypothesized directions. These findings support the FWB-5 Index as a parsimonious measure for examining financial well-being and its associations with health across HRS waves, including those predating the CFPB-FWB Scale.

PhenoAge, a blood-based biological aging clock, is associated with cancer outcomes, a finding taken as evidence that biological aging shapes prognosis. In 3,067 patients undergoing primary tumor resection, PhenoAge acceleration rose with tumor stage. One or two of its eight blood components, differing by cancer type, carried the prognostic signal, and replacing its general-population weights with equal weights did not weaken it. In cancer, PhenoAge may reflect disease burden rather than biological aging.

Purpose: Posterior pharyngeal wall (PPW) thickening is a common structural manifestation of treatment-related internal lymphedema in patients with head and neck cancer and can be evaluated using lateral videofluoroscopic images acquired during modified barium swallow (MBS) studies. However, no standardized method currently exists for quantitative measurement. Methods: This two-phase methodological study first conducted a scoping review to characterize existing approaches to PPW thickness measurement and identify the principal factors contributing to variation in PPW thickness measurement. Eight studies were identified with variation observed across these measurement factors. An empirical reliability study was then performed using 107 MBS studies from 93 Head and Neck Cancer patients. Two independent raters annotated 10 anatomical landmarks on pre-swallow and post-swallow resting frames, and a computer-assisted algorithm quantitatively generated seven candidate PPW thickness measurements based on different vertebral reference points and measurement directions. Reliability was assessed intraclass correlation coefficients (ICCs). Results: Inter-rater reliability was highest for the measurement originating from the midpoint of C3 and extending perpendicular to the anterior C3 surface (ICC, 0.78; 95% CI, 0.72-0.83). Among the normalization references, C2-C4 vertebral length demonstrated greater reliability than C3 vertebral length, with an ICC of 0.88 (95% CI, 0.85-0.91). No significant differences were observed between pre-swallow and post-swallow PPW thickness measurements across the evaluated approaches. Conclusions: These findings provide a methodological foundation for more standardized and reliable quantification of PPW thickness, with potential to expand the role of MBS from the functional evaluation of swallowing to include quantitative characterization of pharyngeal structural abnormalities.

Blood-based neurofilament light (NfL) is associated with all-cause mortality, but whether this reflects neuroaxonal injury or overlaps with cognitive and functional decline remains unclear. We followed 559 individuals aged 92-107 from three oldest-old cohorts for up to 16 years: cognitively healthy Dutch centenarians and two population-based Danish cohorts. Plasma biomarkers (A{beta}42/40, A{beta}40, A{beta}42, pTau-181, pTau-181/A{beta}42, NfL, and GFAP; Dutch cohort; NfL, Danish cohorts), cognition, and activities of daily living (ADL) were evaluated against mortality. NfL was the only biomarker consistently associated with mortality (HR 1.21-1.38), whereas AD-specific biomarkers did not predict mortality despite reflecting post-mortem AD neuropathology. NfL remained associated with mortality independently of cognition and ADL; ADL mediated 13-21% across cohorts and cognition 28%, only where cognitive impairment was widespread. Mortality at extreme age is more closely linked to neuroaxonal injury than to upstream AD pathology, with NfL capturing biological vulnerability beyond cognitive and functional decline.

Objectives: To determine how strongly clinical referral text, rather than image content, drives the decisions of multimodal large language models (LLMs) differentiating glioma from brain abscess on MRI. Materials and methods: This retrospective single-center study included 60 contrast-enhanced T1-weighted brain MRI examinations (30 gliomas, 30 abscesses; reference-standard diagnosis). One axial slice per case was paired with five referral texts (explicitly inflammatory to explicitly tumor-suggestive) under three prompt instructions and submitted 20 times each to GPT-4o, Claude Sonnet 5, Gemini 3.5 Flash, and Qwen 3.6 (72,000 answers). The primary outcome was the text-driven span (difference in tumor-call frequency between explicit inflammatory and explicit tumor text). Two blinded radiologists classified the images without text; sensitivity, specificity, and accuracy were compared descriptively. Results: Patients (29 women) had a mean age of 55.8 +/- 16.4 years. With neutral text, models answered tumor in 45-67% of queries (specificity, 0.38-0.78), versus 52% and 60% for radiologists (specificity, 0.77 and 0.63). The text-driven span was 0.79-1.00 (p < 0.001 for all models), whereas confirmed gliomas and abscesses differed by only 0.01-0.10 in three models and 0.25 in Gemini 3.5 Flash. Subtle wording produced at least 83% of the span in three models. An instruction to ignore the text reduced the span to 0.48-0.70 in three models and to 0.04 in Gemini 3.5 Flash, partly through reversion to a tumor default. Model accuracy was 0.52-0.73 versus 0.78 and 0.73 for radiologists. Conclusion: Multimodal LLMs deferred to referral text over imaging evidence in differentiating glioma from brain abscess, and prompt instructions did not reliably correct this.

Lesion network mapping (LNM) has provided a compelling framework in clinical neuroscience, linking focal brain damage at heterogeneous anatomical locations to shared functional networks across behavioural syndromes. Recently, the methodological specificity of LNM has come under discussion following observations that derived network maps tend to converge onto the normative connectome's generic degree distribution. While current efforts largely treat this convergence as a statistical inference challenge addressable via post-hoc spatial null models or individual-level analyses, we identify a foundational upstream measurement artifact operating within individual patients: the lesion "bystander effect". Because pathological lesions follow biological constraints rather than functional brain parcellations, a typical lesion encapsulates both a symptom-driving critical node and surrounding asymptomatic bystander tissue. By averaging connectivity profiles across the entire lesion mask, LNM mechanically dilutes focal clinical signals. Across empirical clinical cohorts, individual-level machine learning classifiers, synthetic within-cohort lesion unions, and whole-brain spherical simulations, we demonstrate that increasing lesion extent deterministically drives network maps toward the connectome degree vector following a mathematical power law. Fulfilling the clinical promise of LNM will therefore require computational strategies that isolate the causal functional core from surrounding bystander tissue prior to network projection.

Chinese is one of the living languages with the highest homophone densities, making homophone processing a critical component of language comprehension and orthographic decoding, as well as an ideal model for examining its cognitive mechanisms. This study assessed mono- and bi-character homophone selection function in 129 Mandarin- and Cantonese-speaking participants, including 66 cognitively healthy individuals, 17 with semantic variant (sv) primary progressive aphasia (PPA), 23 with the nonfluent/agrammatic variant (nfv) PPA, and 23 with the logopenic variant (lv) PPA. All PPA variants showed impairment in homophone selection, with most significant deficits in svPPA and lvPPA. ROC analyses demonstrated its strong diagnostic utility for distinguishing PPA from controls. Voxel-based morphometry revealed that homophone performance correlated with gray-matter volume in left anterior and inferior temporal regions, with variant-specific patterns observed. These findings described the potential role of homophone selection in diagnosing PPA syndrome in Chinese speakers, illuminate crosslinguistic influences on language impairment in PPA, and provide novel insights into the cognitive processing of homophones.

Background and Objectives Transient ischemic attack (TIA) diagnosis is challenging because it relies mainly on clinical assessment and exclusion of common differential. Functional MRI (fMRI) may provide a diagnostic biomarker by detecting transient brain dysfunction not visible on conventional imaging. This pilot study aimed to characterize acute motor-task BOLD activation in patients with suspected DWI-negative TIA and compare functional patterns between confirmed TIA and alternative diagnoses. Secondary objectives were to assess longitudinal changes in activation and connectivity at 3 months. Methods In this prospective single-center study, adults presenting within 24 h of transient anterior-circulation neurological symptoms underwent conventional MRI and additional motor-task fMRI. Patients with acute DWI lesions were excluded. Final diagnosis was established after 3-month follow-up. Whole-brain analyses, individualized motor cortex region-of-interest (ROI) analyses, and generalized psychophysiological interaction (gPPI) connectivity analyses were performed during the acute phase and at 3 months. Results Eight patients were enrolled; six had complete longitudinal imaging data and were included in the final analyses. Final diagnoses were TIA (n=2), migraine aura (n=2), functional neurological disorder (n=1), and hypoglycemia (n=1). Whole-brain and ROI analyses showed reduced bilateral motor-related BOLD responses during the acute phase in both TIA patients and in the patient with hypoglycemia. In contrast, patients with migraine and functional neurological disorder exhibited less marked longitudinal changes. In TIA patients, motor cortex activation increased at 3 months, tending to normalize toward values observed in non-TIA participants. Connectivity analyses revealed heterogeneous patterns, but the 2 TIA patients demonstrated increased longitudinal connectivity between the motor cortex contralateral to symptoms and other sensorimotor regions, including supplementary motor, premotor, parietal, and postcentral cortices. Discussion These preliminary findings suggest that TIA may be associated with transient bilateral reductions in motor-network activity. Functional abnormalities appeared to improve over time, supporting the hypothesis of reversible network dysfunction. Although limited by the very small sample size, this study provides the first longitudinal assessment of emergency task-based fMRI in rigorously defined TIA and supports further investigation of acute functional MRI as a potential biomarker of transient cerebral dysfunction.

Background: Added sugars during early development may influence adiposity, but whether associations vary by developmental timing and child sex is unclear. Objectives: To determine whether child sex modifies associations of added sugars and fructose exposure from pregnancy through 24 months with body composition at 24 months. Methods: This secondary longitudinal analysis included 130 mother-child dyads. Maternal diet was assessed during mid-pregnancy and child diet at 2wk and 6mo (infancy), and 12mo and 24mo (toddlerhood). Multivariable regression tested sex-by-exposure interactions for anthropometric outcomes; 42 children had usable dual-energy X-ray absorptiometry (DXA) data. Results: Prenatal added sugars were more positively associated with percentage body fat and total and regional fat mass in girls than in boys. Prenatal fructose was more positively associated with waist circumference in boys. During infancy, added sugars were more positively associated with percentage body fat and total and regional fat mass in boys, whereas fructose was more positively associated with fat-free mass in girls. Conclusions: Associations of early-life sugars exposure with body composition varied by developmental timing and child sex. Findings from the smaller DXA subset are exploratory and require replication but suggest that susceptibility to early dietary sugars may differ across developmental windows.

Partitioned-heritability studies have shown which tissues carry the heritability of complex traits, but whether the regulatory layer whose annotations are most informative for heritability within a tissue differs between traits has not been tested, because measured annotations for different layers come from different sets of biosamples. Sequence-to-function models predict the allelic effect of every variant on several regulatory layers with a single model, so that layers can be compared on the same variants within the same tissues. We grouped AlphaGenome Atlas predictions into 12 lineages and six regulatory layers and partitioned the heritability of 11 complex traits with stratified LD score regression. Traits localized largely to the expected lineages, and within these lineages the most informative predicted layer differed between traits (P = 3.6 x 10^-8). Predicted chromatin-accessibility (DNase) effects were the most informative for height and bone mineral density, whereas predicted transcription-initiation (CAGE) effects were the most informative for Alzheimer's disease, most brain traits and cholesterol levels. The difference persisted across sensitivity analyses, and the predicted contrast held for three related traits whose layers had been predicted in advance (P = 2 x 10^-7). Predicted CAGE effects also remained informative after conditioning on fine-mapped eQTLs and measured regulatory annotations. These results identify differences across traits in the heritability information captured by predicted regulatory layers and offer testable hypotheses for investigating the regulatory basis of complex traits.

Background: Commercial spaceflight now carries civilians on short missions, yet their stress-axis physiology is poorly described. Earlier shuttle data suggest that sympathetic responses predominate after short flights and glucocorticoid responses after long ones. The neutrophil-to-lymphocyte ratio (NLR) is a low-cost readout of both. Methods: Preregistered secondary analysis of open clinical laboratory data (NASA Open Science Data Repository OSD-569 and OSD-575) from the four Inspiration4 crew members, each sampled seven times: three draws before launch, one the day after return (R+1) and three during recovery (R+45 to R+194). Log-transformed NLR was modelled with a linear mixed model with a random intercept per crew member and Kenward-Roger inference. Lymphocyte and eosinophil counts and C-reactive protein were secondary outcomes with false discovery rate correction. Results: All 28 draws were analysed. NLR at R+1 did not differ from pre-flight (fold change 1.04; 95% CI 0.81 to 1.34; p = 0.73); two crew members increased (+25% and +46%) and two decreased (-27% and -11%). During recovery NLR was 23% lower than pre-flight (fold change 0.77; 95% CI 0.64 to 0.92; p = 0.005). No secondary outcome changed after correction; lymphocytes were 18% lower at R+1 (95% CI -35% to +4%; q = 0.31). Between-person differences explained 58% of the variance. An ordinary regression ignoring repeated measures missed the recovery effect (p = 0.08). Conclusions: A 3-day orbital flight produced no consistent leukocyte stress signature one day after return. The lower NLR during recovery suggests that pre-flight values, not recovery values, deviate from the resting state, which has implications for how baselines are defined in short-duration spaceflight studies.